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LONG-TERM  GOALS 

My  project  addresses  statistical  and  stochastic  problems  in  the  following  fields:  Lagrangian  data 
assimilation  (1),  Lagrangian  prediction  (2),  development  and  validation  of  Lagrangian  stochastic 
models  (LSM)  (3).  The  long  range  scientific  objectives  of  this  study  comprise  rigorous  determining 
limits  of  predictability  for  Lagrangian  motion  in  semi-enclosed  seas,  littoral  zones,  and  straits  on  time 
scales  of  days  and  weeks,  elaborating  concrete  prediction  schemes,  developing  optimal  Lagrangian 
data  assimilation  algorithms,  and  identification  of  multi  particle  stochastic  models  aimed  at 
incorporating  them  to  ocean  circulation  models  (OCM). 

OBJECTIVES 

The  objectives  for  the  third  year  of  research  were: 

•  Development  and  verification  of  estimation  algorithms  for  Eulerian  velocity  field  given 
Lagrangian  data  and  distributed  tracer  observations  based  on  the  non-linear  filtering  approach. 

•  Construction  and  verification  of  data  fusion  procedures  for  computing  surface  velocities  given 
a  circulation  model  output,  drifter  data  and  observations  of  a  continuously  distributed  tracer. 

•  Extension  of  the  Lagrangian  subgridscale  model  [1]  to  multi  particle  motion  for  improvement 
of  Lagrangian  transport  prediction. 

•  Investigation  of  particle-pair  statistics  in  the  framework  of  a  LSM  with  autocorrelated  forcing. 

APPROACH 

I  develop  theoretical  approaches  to  the  Lagrangian  prediction,  Lagrangian  data  assimilation  as  well  as 
to  data  fusion,  in  context  of  the  theory  of  random  processes  and  fields  covered  by  stochastic  partial 
differential  equations.  I  design  computational  algorithms  derived  from  the  theoretical  findings.  During 
the  last  year  I  also  made  use  of  possibility  theory  and  fuzzy  logic  to  design  procedures  for  optimal 
combining  (fusion)  data  coming  from  different  sources.  A  significant  part  of  validating  the  algorithms 
is  a  testing  them  via  stochastic  simulations.  Such  an  approach  provides  us  with  an  accurate  error 
analysis.  Together  with  my  collaborators  from  Rosenstiel  School  of  Marine  and  Atmospheric  Research 
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WORK  COMPLETED 

1.  Assimilation  and  data  fusion. 

The  problem  of  assimilation  both  drifter  and  image  data  in  an  Eulerian  model  is  formulated  as  an 
optimal  filtering  problem.  A  quasi-optimal  solution  of  the  problem  has  been  found  for  a  wide  class  of 
stochastic  models  which  are  able  to  realistically  mimic  the  upper  ocean  turbulence.  A  theoretical  error 
analysis  is  developed  for  some  important  asymptotical  cases.  The  approach  requires  a  particular 
knowledge  of  the  Eulerian  velocity  field  statistics.  For  this  purpose  special  methods  were  developed  to 
estimate  appropriate  parameters  from  Lagrangian  data. 

For  the  same  problem  a  different  approach  was  developed  and  tested  based  on  the  possibility  theory 
ideas,  [2,3]. 

2.  Lagrangian  prediction. 

A  procedure  was  developed  to  reduce  errors  in  the  relative  dispersion  obtained  from  a  numerical  model 
by  employing  drifter  statistics  estimated  from  observational  data  sets.  The  method  is  based  on  the  PI 
results  regarding  to  the  relative  dispersion  in  multi  particle  Lagrangian  stochastic  models  [4],  Earlier  a 
similar  approach  resulted  in  improvement  of  modeling  the  absolute  dispersion  [1]. 

3.  Development  of  multi  particle  Lagrangian  stochastic  models. 

A  new  stochastic  model  for  describing  a  particle  pair  has  been  formulated  where  the  acceleration  has  a 
finite  correlation  time.  Such  an  approach  allows  for  consideration  of  the  particle  pair  statistics  as 
functions  of  both  Kubo  and  Stokes  numbers.  A  comprehensive  analytical  and  Monte  Carlo 
investigation  of  the  Lyapunov  exponent  was  carried  out  to  establish  the  dependence  of  predictability 
and  mixing  properties  on  these  two  parameters. 

RESULTS 


1 .  The  first  principal  result  in  the  assimilation  problem  is  derivation  and  verification  of  optimal 
filtering  equations  for  the  underlying  Eulerian  velocity  given  a  model  output  u m(t,  x),  observations  of 
Lagrangian  trajectories  rft)\,  k  =  1,2,...  driven  by  the  real  velocity  field,  and  grid  observations  of  a 
continuously  distributed  tracer  c(t,x)  satisfying  a  standard  advection-diffusion  equation  with  a  white 
noise  forcing.  For  simplicity  we  show  the  filtering  equations  in  the  case  of  one  drifter  and  tracer 
observations  in  a  single  point  x  only 


du  =  (— u/t  +  F(i,  x))dt  +  -B(x,  r(t))B(r(t),r(t))  1(dv  +  (v/r  —  F(i,  r (t))dt) 
— rVcR(x,  x)_1(dc  +  (u  •  Vc  —  DV2c)dt ) 


(1) 


dT/dt  =  —  2Y/t  +  .B(x,  x)  —  £>(x,  r(t))B(r(t),  r(t))  1B(x,  r(t))T  —  TVcR(x,  x)  1VcTr7  (2) 
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Here  u  is  the  velocity  estimate,  F  its  error  matrix,  F(/,  x)  =  dumldt  -  u Jx,  x  the  Lagrangian 
correlation  time,  D  diffusivity,  R(x,y)  space  covariance  of  the  forcing  and  B(xy)  the  space  covariance 
of  velocity  fluctuations.  In  addition  we  have  developed  a  procedure  of  estimating  the  crucial 
parameters  of  the  covariance  B(xy)  from  Lagrangian  data. 

Essentially,  (1-2)  is  a  Kalman  filter  written  in  the  continuous  form  which  has  been  applied  before  for 
data  assimilation  (e.g.  [5,6]).  The  news  here  is  that  for  a  formulated  general  enough  model  it  gives  an 
exact  solution  rather  than  an  approximation  based  on  a  linearized  version  of  a  sophisticated 
hydrodynamics  model.  We  have  tested  the  discretized  version  of  this  algorithm  modified  for  multiple 
drifters  in  the  case  with  no  tracer  observations.  Comparison  of  the  error  with  the  simplest  interpolation 
algorithm  (SA),  [7],  for  different  values  of  the  space  correlation  radius,  shows  that  the  perfonnance  of 
the  new  method  is  clearly  better,  Fig.  1 


Figure  1.  Left  panel:  velocity  estimation  error  for  algorithm  (1,2)  vs  time  (days)  at 
different  values  of  space  correlation  radius  R=2  km  (square),  4  km  (cross),  6  km  (star), 

8  km  (circle),  10  km  (plus).  Right  panel:  same  for  SA} 

We  also  have  tested  another  approach  to  the  same  problem  based  on  the  possibility  theory  [1,2].  In  the 
simplest  case  of  an  isotropic  confidence  region  for  a  circulation  model  output  and  tracer  observations 
only  ,  the  fusion  formulas  are  as  follows 

^  A(Aum  +  Bvm  +  E )  ^  B(Aum  +  Bvm  +  E )  .  . 

u  =  u m - A*  +  B2 - ’  V  =  Vm - ATTb~2 -  (3) 

where  um,  vm  are  the  velocity  estimates  obtained  from  the  circulation  model  (possibly  with  drifter  data 
assimilation)  and  A  =  dc/dx,  B  =  dc/dy,  E  =  dc/dt  are  taken  at  the  underlined  point. 

This  method  allows  a  rigorous  error  analysis.  One  of  our  findings  is  that  if  <Tmodei  =  (|u,„  -  u|2)  is 
small  enough  where  u  is  the  real  velocity  and  the  angle  brackets  mean  the  space  averaging,  then  cr fusion 
(|  um  -  u|2)  is  twice  as  less  as  <Tmodei- 

To  test  the  algorithm  we  used  an  idealized  velocity  field  in  a  rectangle  region  which  is  a  superposition 
of  the  gyre  with  the  stream  function  p  =  co(x  +y  )  and  periodic  ageostrophic  perturbations.  It  is 
assumed  that  some  particular  model  is  not  able  to  capture  those  perturbations  and  its  output  coincides 
with  p .  Next  a  tracer  is  spread  over  the  region  by  the  real  velocity.  The  estimation  result  for  the  zonal 
component  after  one  day  is  shown  in  Fig.  2.  The  resulting  crfusion  =  7%  is  in  full  agreement  with  the 
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theoretical  error  analysis  since  we  have  taken  (fniode,  =  10%,  but  what  is  more  important  that  the  method 
is  capable  to  capture  the  qualitative  structure  of  the  current. 


Figure  2.  Real  u  (left),  model  u  (center)  and  the  estimated  velocity  by  fusion  of  model 
output  and  tracer  observations  via  (3)  (right) 


2.  The  most  important  finding  in  our  Lagrangian  prediction  study  is  a  computationally  efficient 
algorithm  of  correcting  model  generated  trajectories  in  a  way  such  that  the  resulting  relative  dispersion 
is  close  to  that  of  estimated  from  observations.  For  simplicity  we  show  only  the  corresponding 
equation  for  the  zonal  component  of  Lagrangian  velocity.  Namely,  let  uft,  r)  and  um(t,  r)  be  a  real  and 
model  velocities  described  for  any  label  r  by  the  Langevin  equation  with  different  parameters, 

(< <J„  T„)  and  (<Jm,  Tm)  respectively.  We  have  proven  that  under  certain  realistic  conditions  the  corrected 
velocity  u  =  um+Tj  generates  the  same  absolute  and  relative  dispersion  if  the  missing  component  77 
satisfies  the  following  stochastic  differential  equation 


drj  dum 

—  =  0,— — 
dt  dt 


+  bum  +  cry 


where  the  coefficients  a,  b,  c  are  expressed  in  terms  of  (<7r,  T,)  and  (am,  Tm)  exactly  as  in  [1]  and  um  = 
Lum  where  A  is  a  linear  operator  applied  to  the  label.  L  depends  on  some  parameters  characterizing  the 
space  covariance  structure  of  the  real  and  model  velocity  fields.  In  particular  if  the  parameters  are  close 
for  the  real  and  model  velocity  fields,  then  L  is  the  identical  operator  and  the  method  reduces  to  the 
known  one  [1],  where  only  one-particle  statistics  are  intended  to  match.  The  mentioned  parameters 
such  as  the  space  correlation  scale  and  the  energy  slope  on  large  scales  are  efficiently  estimated  from 
Lagrangian  data  [8]. 


3.  In  the  new  suggested  stochastic  particle -pair  model,  the  dependence  on  Stokes  (£j  and  Kubo  (h) 
numbers  of  the  Lyapunov  exponent  (LE)  of  particle  trajectories  reveals  the  presence  of  a  region  in 
parameter  space  ( e. ,  h)  where  LE  changes  sign,  thus  signaling  the  transition,  Fig.  3. 
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Figure  3.  Left  panel:  dependence  of  LE  on  Kabo  number  (h)  for  different  Stokes 
numbers  (e).  Right  panel:  dependence  of  LE  on  Stokes  numbers  (e)  for  different  Kubo 

number  (It) 

In  particular,  for  large  Stokes  number,  a  regime  characterized  by  the  fonnation  of  shocks,  we  found  a 
chaotic  region  in  parameter  space  (e,  h ),  where  LE  becomes  positive.  Inertia  is  therefore  responsible 
for  a  transition  from  a  strong  clustering  regime,  originated  by  the  compressible  nature  of  the  flow,  to  a 
chaotic  regime.  The  latter  is  observed  in  a  range  of  h  such  that  the  time  correlation  of  fluid  gradients  is 
long  enough  to  provide  substantial  stretching,  but  not  too  long  to  get  particles  trapped  in  compressing 
regions. 

IMPACT/APPLICATIONS 

1 .  The  suggested  data  fusion  procedures  for  estimating  Eulerian  velocity  given  Lagrangian 
observations  ,  images,  and  a  model  output  could  be  a  supplement  to  existing  assimilation  algorithms 
explicitly  involving  OGCM. 

2.  The  suggested  procedure  of  adding  a  missing  component  to  Lagrangian  velocities  coming  from 
OGCM  can  potentially  improve  their  perfonnance  for  particle  pair  statistics  such  as  the  relative 
dispersion,  Lyapunov  exponent  and  finite  time  Lyapunov  exponent. 

3.  The  developed  theory  of  the  Lyapunov  exponent  along  with  estimates  of  the  Stokes  and  Kubo 
numbers  for  the  oceanic  conditions  would  lead  to  classification  different  regions  with  respect  to  the 
mixing  rate  and  the  accuracy  of  prediction  of  the  Lagrangian  trajectories. 

TRANSITIONS 

The  developed  trajectory  correction  algorithm  was  used  in  RSMAS  to  test  it  in  NCOM  circulation 
model.  Its  new  particle  pair  version  is  planned  to  implement  in  an  OGCM  by  the  same  RSMAS  group. 
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